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In this chapter we provide an overview of the current status of the simulations and modelling of 
the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated 
by the characteristic scales of the problem and the SKA instrumental properties and the planned 
survey parameters. Current simulations include most of the relevant physical processes. They can 
follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, 
although small-scale physics still needs to be included at sub-grid level. However, despite a signif¬ 
icant progress in developing novel numerical methods for efficient utilization of current hardware 
they remain quite computationally expensive. In response, a number of alternative approaches, 
particularly semi-analytical/semi-numerical methods, have been developed. While necessarily 
more approximate, if appropriately constructed and calibrated on simulations they could be used 
to quickly explore the vast parameter space available. Further work is still required on including 
some physical processes in both simulations and semi-analytical modelling. This hybrid approach 
of fast, approximate modelling calibrated on numerical simulations can then be used to construct 
large libraries of reionization models for reliable interpretation of the observational data. 
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1. Introduction 

The Cosmic Dark Ages, Cosmic Dawn and Epoch of Reionization, span altogether about 
a billion years between the CMB last scattering surface at redshift z ~ 1100 and the complete 
ionisation of hydrogen in the inter-galactic medium (IGM) at z ~ 6. This phase is quite distant 
from us and difficult to study, remaining one of the last poorly understood epochs of the Universe. 
During this time, the very first stars and galaxies formed, then gradually ionized the intergalactic 
medium and enriched it with metals, thereby laying the foundations for the cosmic structures we 
see today. Better knowledge of these epochs is therefore of key importance for understanding 
galaxy formation and evolution. 

The main obstacle to further progress is the scarcity of observational data, which currently 
mostly probes the tail-end of reionization (both Ly-a source surveys and IGM absorption lines 
probe low neutral fractions, due to the high optical depth to resonant scattering of such radiation by 
even small amounts of neutral hydrogen) or are integral measures of its history (Cosmic Microwave 
Background, CMB, optical depth and polarization; kinetic Sunyaev-ZeTdovich effect, kSZ; Near 
Infrared Background, NIRB). The redshifted 21-cm signal promises to provide full 3D tomographic 
observations of the intergalactic medium throughout, and possibly even before reionization. The 
first generation of experiments, currently ongoing, will most likely provide just a statistical detec¬ 
tion of the signal, e.g. power spectra, variance and/or probability distribution functions (PDFs) 
and their higher moments. In contrast, SKA should be able to also perform imaging and 3D map¬ 
ping, as well as go much deeper in redshift, due to its far superior sensitivity, thereby completely 
transforming reionization research. 

The shortage of observational constraints has meant that simulations have played and continue 
to play a larger role than in other areas, since we need to rely on modelling to understand even 
the basic properties of the expected signals. This in turn steers the design of the observational 
experiments targeted at this science. Furthermore, better understanding of the expected 21-cm 
signatures and their cross-correlations with other observations is important for reliably confirming 
any detection, since the signal is weak compared to the foregrounds which are several orders of 
magnitude stronger. Separating the signal and the foregrounds will not be trivial (see e.g. Chapman 
et al. 2015, reference: PoS(AASKA 14)005) and subtraction may leave contaminating residuals 
mimicking the signal. The simulations will have a further important role to play in the interpretation 
of any detections. The reason for this is the disconnect in scales between the objects we want to 
study, i.e. the first galaxies, which are very small and faint, thus largely impossible to observe 
directly, and the actual 21-cm and other signals we will detect, which typically relate to large-scale 
patterns such as the reionization patchiness. The latter arc caused by large number of clustered, 
individually dim sources. Therefore, the interpretation of any signals in terms of the nature and 
properties of the first galaxies is non-trivial and requires detailed modelling. 

In recent years two basic approaches emerged in modelling the reionization patchiness - ei¬ 
ther by direct numerical simulations (e.g. Gnedin & Ostriker, 1997; Ciardi et al., 2000; Gnedin, 
2000a; Ricotti, Gnedin & Shull, 2002; Iliev et al., 2006, 2014) or through semi-analytical/semi- 
numerical modelling (e.g. Furlanetto, Zaldarriaga & Hernquist, 2004; Zahn et al., 2007; Alvarez 
et al., 2009; Mesinger & Furlanetto, 2007a; Choudhury, Haehnelt & Regan, 2009; Santos et al., 
2010; Battaglia et al., 2013), along with some intermediate methods involving simplified simula- 
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Figure 1: Spatial slices from radiative transfer simulations of reionization. Shown are the ionized (or¬ 
ange/blue) and neutral gas density (green/pink) at box-averaged ionized fractions by mass of x m ~ 0.4 for 
(left) a small, high resolution volume (47 Mpc/h, 612 3 radiative transfer grid) and a very large volume, com¬ 
parable to the expected full FOV for SKA1-LOW (500 Mpc/h, corresponding to ~ 5° x 5° on the sky, 250 3 
grid) (Dixon et al., in prep.). 

tions (Kohler, Gnedin & Hamilton, 2007; Thomas & et al., 2009). These different lines of attack are 
in practice highly complimentary and each has its advantages as well as disadvantages, as discussed 
in more detail in the next section. 

2. Overview of current modelling 

2.1 Radiative transfer simulations 

The simulations of cosmic reionization fall in two broad types, as a consequence of the huge 
range of scales involved, which in practice cannot be covered in a single simulation on current 
computer hardware. The first type consists of small volume, high-resolution simulations which 
can be used to study in detail the formation of early galaxies and their radiative and supernova 
feedback on the gas with fully-coupled radiative hydrodynamic simulations. Such simulations can¬ 
not capture a large enough volume to represent the global reionization process or its observational 
signatures (e.g. Gnedin & Ostriker, 1997; Ciardi et al., 2000; Gnedin, 2000a; Ricotti, Gnedin & 
Shull, 2002). The second type consists of large-scale simulations, which instead follow volumes 
sufficiently large to study the global evolution, but lack resolution to directly resolve the small- 
scale physics (e.g. Iliev et al., 2006; Trac & Cen, 2007; Iliev et al., 2014) (see Figs, [l] and |j). This 
type of simulations decouple the radiative transfer, star formation and feedback processes from the 
underlying structure formation. The latter is done first, typically using large-scale N-body sim¬ 
ulations, which is then followed up by detailed radiative transfer and non-equilibrium chemistry 
simulations. This approach allows for much larger dynamic range, and thus the larger volumes 
required for understanding the global reionization process and its observational signatures. In this 
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Figure 2: Space-frequency slice of the 21-cm emission from the intergalactic medium showing the evolution 
of the reionization geometry over time, from fully-neutral (blue) to mostly ionized (green). This is based 
in a radiative transfer simulation with 500Mpc/h volume and 300 3 cells. The color shows the differential 
brightness temperature in mK at the full simulation resolution. (Dixon et al., in prep.) 


case the detailed physical processes at unresolved scales have to be included using subgrid recipes, 
which are themselves based on high-resolution simulations or other modelling (e.g. Iliev et ah, 
2007; McQuinn et al., 2007; Ahn et ah, 2012, 2014). Both types of simulations have played a 
major role in the significant recent advances in this area. We now have a good understanding of the 
characteristic scales of the process (Friedrich et ah, 2011; Iliev et ah, 2014), the physical processes 
affecting the patchiness and the various observational features which might help to discriminate 
between different EoR models (e.g. Wise & Abel, 2008; Baek et ah, 2009; Ahn et ah, 2012; Iliev 
et ah, 2012; Alvarez & Abel, 2012; O’Leary & McQuinn, 2012; Park et ah, 2013; Jensen et ah, 
2013; Sobacchi & Mesinger, 2014). 

2.2 Semi-numerical simulations 

Some of the early approaches to characterise the 21cm signal from the epoch of reionization 
have relied on more simple analytical models, mostly based on the excursion-set approach. These 
can be very useful to quickly generate the signal, in particular the power spectrum of 21-cm bright¬ 
ness temperature fluctuations (Furlanetto, Zaldarriaga & Hernquist, 2004; Furlanetto, McQuinn & 
Hernquist, 2006; Sethi, 2005). The speed and ease of use, allows for straightforward tests of the 
ability of a given experiment to constraint cosmological and astrophysical parameters (Santos & 
Cooray, 2006; McQuinn et ah, 2006; Mao et ah, 2008) and more generally to rapidly explore the 
vast parameter space available. The analytical models have also been useful to understand the pos¬ 
sible contributions to the 21 cm signal at high redshifts (Barkana & Loeb, 2005; Pritchard, 2007), 
although it is at low redshifts (z < 10) that they seem to provide a better description of the 21-cm 
2-point correlation function (Santos et ah, 2008). However, these models have several issues in 
properly dealing with the spatial distribution of the reionization process, such as bubble overlap 
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and ignore complicated astrophysics during reionization. Such complexities are better handled in 
full (but expensive) numerical simulations, as was discussed above. 

An intermediate approach that has become more popular are the so called semi-numerical 
21cm simulations which try to merge the speed and ease of use of analytical models with the ability 
to follow the spatial evolution of reionization in more detail as provided by numerical simulations. 
The basic feature of these semi-numerical simulations is that they start from a random Gaussian 
field realization of the cosmological density field, or a full N-body simulation data, but then replace 
the time-consuming radiative transfer (RT) with an excursion-set approach (Furlanetto, Zaldarriaga 
& Hernquist, 2004) that determines if a region is ionized by comparing the number of ionizing 
photons to the number of atoms (plus recombinations) inside it. As in conventional cosmological 
RT simulations, this algorithm can be applied to discrete halo source fields, obtained either from 
traditional N-body simulations (e.g. Zahn et al., 2007; Choudhury, Haehnelt & Regan, 2009; Zahn 
et ah, 2011) or through faster methods involving excursion-set formalism and perturbation theory 
(Mesinger & Furlanetto, 2007a; Santos et ah, 2010; Mesinger, Furlanetto & Cen, 2011; Alvarez 
et ah, 2009; Geil & Wyithe, 2008). The latter approach results in some reduction in accuracy 
(see Fig. j|; Zahn et ah, 2011; Majumdar et ah, 2014), but substantially increases the achievable 
dynamic range. Some of these approaches have also been extended to simulating the earlier epochs 
of the 21cm signal: the Cosmic Dawn/Dark Ages (Mesinger, Furlanetto & Cen, 2011; Santos et ah, 
2010). During the Dark Ages, the 21cm signal is governed by soft UV and X-ray photons, which 
can have mean free paths of hundreds of Mpc, requiring large volumes, which makes it a challenge 
to implement this in a complete radiative transfer/gas dynamics algorithm. 

Semi-numerical simulations of reionization have been extensively tested against numerical 
simulations, with good agreement on moderate to large scales (>Mpc), relevant for most SKA sci¬ 
ence (Mesinger, Furlanetto & Cen, 2011; Zahn et ah, 2011; Majumdar et ah, 2014). The approxi¬ 
mations start to break down however on non-linear scales, and care should be taken by comparing 
against more accurate methods when extending their application to any non-standard problems. 
Nevertheless, their speed and efficiency allows for rapid exploration of the large-parameter space 
of uncertainties, ushering in an exciting era of 21cm astrophysical parameter studies (e.g. Mesinger, 
Ewall-Wice & Hewitt, 2014; Pober et ah, 2013). 

The analysis of Zahn et ah (2011) has been largely focused on the comparison of the spheri¬ 
cally average power spectrum of the 21-cm signal from different simulations. Their study suggested 
that these semi-numerical techniques can mimic the 21-cm power spectrum from that of a radiative 
transfer simulation with > 70% accuracy at length scales relevant for the present and future sur¬ 
veys. This study was extended in more recent works (Mesinger, Furlanetto & Cen, 2011; Majumdar 
et ah, 2014) to include the comparison of EoR history, the redshift space anisotropies in the power 
spectrum and the morphology of the ionization fields (Fig. |3j). Those analyses suggest that the 
halo based semi-numerical simulations can predict the history of EoR with an accuracy of > 90% 
when compared with radiative transfer simulations for length scales k < 1.0Mpc” 1 , whereas the 
same accuracy for the Press-Schechter based technique is > 40% at same length scale range. They 
have also shown that the anisotropies in the power spectrum of the 21-cm signal due to redshift 
space distortions can also be mimicked by both type of the semi-numerical simulations with an 
accuracy of 85% or more, which is well within the noise uncertainties of ~ 2000 hi - LOFAR or 
~ 1000 hr SKA1-LOW observation at 150 MHz, provided that the actual peculiar velocities are 


5 



EoR modelling for SKA 


Ilian T. Iliev 



Mpc 


Mpc 


Mpc 


hydro + N-body + RT 


DexM (with halos) 





21cmFAST (without halos) 


Figure 3: Reionization morphology through brightness temperature maps in redshift space for radiative 
transfer simulation (left) vs. two semi-numerical approaches - one which takes the halo positions into ac¬ 
count (middle) and one that does not (right). For details see the original papers Majumdar et al. (2014) (top) 
and Mesinger, Furlanetto & Cen (2011) (bottom). 


used to introduce the redshift space distortions rather than using a perturbative technique. It is also 
found that the cross-correlation between the ionization fields from the halo based semi-numerical 
technique and the radiative transfer simulations is more than 85% during almost the entire period 
of EoR, for length scales k < 1.0Mpc -1 , whereas the same correlation between the conditional 
Press-Schechter formalism and the radiative transfer simulations, while good in some regimes, can 
reach values as low as ~ 10%, specifically at late stages of EoR for the same length scale range. 

These comparisons however, arc also limited as the reionization scenario/model that has been 
compared in these cases is relatively simplistic in nature. The reionization history and the nature 
of the 21-cm signal from this epoch may differ depending on various other factors. One such 
important factor is the effect of radiative feedback on the star formation in low mass halos (< 
1O 9 M 0 in mass). This will most likely affect their star formation efficiency, although the details 
remain unclear (Couchman & Rees, 1986; Rees, 1986; Efstathiou, 1992; Thoul & Weinberg, 1995, 
1996; Gnedin, 2000b; Kitayama et al., 2000; Dijkstra et al., 2004; Hoeft et al., 2006; Okamoto, 
Gao & Theuns, 2008). There has been some effort to include this in numerical (Iliev et al., 2007, 
2012; Ahn et al., 2014) as well as in semi-numerical simulations (Sobacchi & Mesinger, 2013). 
Similarly, another such important issue is the effect of enhanced recombination in the Lyman Limit 
Systems and other small-scale structures, which are not resolvable in any of these simulations and 
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thus require detailed sub-grid modelling (Ciardi et al. 2006; Choudhury, Haehnelt & Regan 2009; 
Sobacchi & Mesinger 2014, Koda et al. in prep., Shukla et al., in prep.). Also, the effect of X-ray 
heating on the spin temperature evolution at the very early stages of reionization (Mesinger, Ferrara 
& Spiegel, 2013; Ghara, Choudhury & Datta, 2014) is another important factor which has not been 
included in these comparison studies. 

3. SKA simulations 

3.1 Basic Simulation Requirements 

The volume and resolution required for proper numerical modelling of the reionization pro¬ 
cess are dictated by both its intrinsic characteristic scales (especially the typical sizes of ionized and 
neutral patches) and the instrumental parameters (FOV, beam size and bandwidth). The simulation 
volume should be large enough to faithfully sample all relevant scales, while the numerical resolu¬ 
tion should be such that every beam/bandwidth is sampled with a sufficient number of elements so 
as to avoid discretization effects and other numerical artifacts. 

The nominal SKA 1-LOW Epoch of Reionization survey at the frequencies relevant for reion¬ 
ization is proposed to have FOV (FWHM) of ~ 3 (at v = 200 MHz) to 10 degrees (at v = 50 MHz), 
which corresponds to ~ 500 Mpc-1.5 Gpc, with maximum resolution of ~ 1 arc min. or roughly 
0.5 Mpc. Therefore simulations of the full SKA1-LOW FOV with sufficient resolution should have 
volume of at least 500 Mpc and grid sizes of at least several thousand cells per side. These simu¬ 
lation parameters arc now just becoming achievable on current hardware. These arc of course only 
very rough estimates and they only relate to the 21-cm sky maps and statistical studies (e.g. power 
spectra, or PDFs). For other purposes, for example 21-cm absorption studies against bright radio 
galaxies (see e.g. Ciardi et al. (2015) reference: PoS(AASKA14)006), a much higher simulation 
resolution will be required, corresponding to frequency channels of ~ kHz, or cells of ~ 10 kpc 
comoving. This is unrealistic when combined with ~ 500 Mpc volumes unless adaptive-mesh re¬ 
finement (AMR) techniques are employed, or else smaller simulation volumes have to be used for 
such work. 

Apart from the purely instrumental requirements discussed above, the simulation parameters 
arc also dictated by the characteristic scales of the reionization itself and the various processes that 
occur at different stages and will be of interest to study. Before the first UV sources form, the 21-cm 
fluctuations arc dictated by the underlying density fluctuations, the temperature of the IGM, as well 
as additional effects like baryon-dark matter displacement (Tseliakhovich & Hirata, 2010, see also 
Maio et al. 2015, reference: PoS(AASKA14)009) and star formation suppression in minihaloes due 
to Lyman-Werner bands photons. Once the first ionizing sources appeal - , they propagate ionization 
fronts into the IGM initially forming small, local HII regions, which then quickly percolate and 
merge locally into larger ones. The same first sources also produce copious amounts of soft-UV ra¬ 
diation and likely some X-rays, which heat the IGM gas and decouple its spin temperature from the 
CMB, making it detectable. Outside dense regions/very high redshifts, the 21-cm line decoupling 
from the CMB also requires sufficiently strong Lyman-a background. All these different processes 
impose some characteristic scales on the 21-cm fluctuations. These are typically long-range mod¬ 
ulations ranging from tens of Mpc (Lyman-a), to ~ 100 Mpc (baryon-dark matter displacement 
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Figure 4: (left) How common are halos of different mass at high redshift? The lines show the mass vs 
redshift for v — a fluctuations in the Gaussian statistics, where v = 8 c /o(z = 0,M) = 1 (i.e. M*, typical 
halo; black, solid), v = 2 (red, dotted), v = 3 (green, short-dashed) and v = 5 (blue, long-dashed) (Iliev 
et al., 2012). (right) The halo bias, f?jj h = Ahh/A p , at redshift z = 9 for three cosmological volumes at 
different resolution, as labeled. Lines are for haloes binned by decades of mass (bottom to top curve) 
10 5 M 0 <M hal0 < 1O 6 M 0 , 10 6 M e <M hal0 < 1O 7 M 0 , 1O 7 M 0 <M halo < 1O 8 M 0 , 10 8 M 0 <M halo < 1O 9 M 0 , 
1O 9 M 0 <M ha i 0 < 1O 1O M 0 , 1O 1O M 0 < Mh a i 0 < 1O U M 0 , 10 n M 0 <M halo < 1O 12 M 0 , and 1O 12 M 0 < M halo . 


and Lyman-Werner) and up to hundreds of Mpc (X-rays). Later-on, likely below redshift ~ 10 
when the reionization becomes more widespread and is driven by larger, atomically-cooling halos, 
the 21-cm fluctuations due to these early backgrounds become less significant, baryon-dark matter 
displacement decreases and the dominant process determining the 21-cm fluctuations become the 
ionization patchiness. Its typical scales and geometry depend on the abundance and clustering of 
the main sources driving the process. Low-mass galaxies are exponentially more abundant, but are 
also liable to get their star formation suppressed by the Jeans mass filtering due to photoheathing, 
as well as other radiative and mechanical feedback mechanisms. The details of these processes are 
still not firmly established and are currently subject of active research. 

The ionized regions growth continues throughout reionization (see Fig. [2] for an example). The 
characteristic patch sizes are dictated by the abundance, clustering and typical luminosities of the 
dominant ionizing sources, as well as various effects which modify the growth like recombinations 
in the IGM gas, Lyman-limit systems and other photon sinks. The intrinsic source clustering results 
from the statistics of the initial Gaussian random noise density fluctuations and depends strongly on 
redshift and the typical mass of the halos hosting the dominant sources. In ACDM the cosmological 
structures form hierarchically, with the smallest ones forming first and then growing and merging 
over time to form larger ones. Consequently, low-mass galaxies likely dominated the ionizing 
photon output throughout reionization. Figure [4] (left) shows the relative abundance of halos of 
different mass vs. redshift. At very high redshifts all halos are rare (v>l) and only below redshift 
z ~ 15 the low-mass galaxies (M = 10 8 — 1O 1O M 0 ) become somewhat more common (v = 2 — 3). 
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Figure 5: (left panel) Probability distribution function per logarithmic radial bin, Rjf r spherical regions 
with radius R as given by the Spherical Average method (Zahn et al., 2007) based on the ionized distribution 
given by a reionization simulation with volume of 607 Mpc and 504 3 grid (Iliev et al., 2014). Shown are 
results for different stages of the reionization process, for ionized fraction by mass x m = 0.1 to 0.8 (left to 
right), (right panel) Dimensionless power spectra, A x , of the volume-weighted ionized fraction, for H II 
regions extracted from the same reionization simulation and at the same EoR stages (bottom to top). 

Consequently, these halos are strongly clustered, with bias with respect to the underlying mass 
distribution well above 1 (Figure ||, right). This strong clustering yields quick percolation of the in¬ 
dividually small Fill regions around each relatively weak source and thus the rapid growth to much 
larger scales. Furthermore, there is a notable power in the density fluctuations at fairly large scales, 
which required large simulation volumes to account for properly (Mesinger & Furlanetto, 2007b; 
Iliev et al., 2014). This additional power means that local photon output is modulated significantly, 
which is not modelled correctly in volumes smaller than about ~ lOOMpc/h per side, resulting in 
artificial suppression of the redshifted 21-cm fluctuations, among others. Different measures of the 
characteristic scales of the patchiness all yield typical peak sizes of tens of comoving Mpc (see Fig¬ 
ure [|), which corresponds to tens of arcmin (see also Sobacchi & Mesinger, 2014). In summary, 
accurate modelling of all these physical processes again points to requirements of large, at least 
several hundred Mpc per side, volumes and sufficient resolution to reliably identify, either directly 
or sub-grid, the main sources and sinks driving the reionization process. 

3.2 Computational Requirements 

As mentioned above, reliable guidance and interpretation of the observations in terms of de¬ 
riving the properties of the first galaxies requires the creation of large libraries of models sampling 
the available parameters space. At present a few full radiative transfer codes can be scaled up to the 
volumes and resolutions required for full SKA1 EoR simulation discussed in the previous section. 
This can only be achieved using massively-parallel approach (in some cases using accelerators 
like Graphical Processing Units, GPUs, and Intel Phi multicore processors) on the largest avail- 
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able computers and a single simulation of this size requires up to tens of millions of core-hours. 
Even with the expected advances in computing technology it is unlikely that a large number of 
such simulations could be practically performed. It is therefore most likely that more approximate 
methods like semi-numerical modelling would have to be employed, guided by fewer detailed full 
simulations. 

Both approaches also require large amounts of data storage. The large N-body structure forma¬ 
tion simulations (needed for precision and accuracy even in semi-numerical models, as explained 
in § |2.2| ), alone require up to several PB of storage per simulation, with further significant storage 
required for the reionization data itself. Accessing and using this data within the community will 
require also investment in databases and other efficient methods for data sharing. 

4. Summary 

We presented an overview of the current status and future prospects for simulations and mod¬ 
elling of the Epoch of Reionization, with focus on the specific requirements for SKA 1-LOW. As 
discussed in detail, extensive modeling is particularly important for Cosmic Dawn and EoR science 
compared to other areas of study because as yet there is very limited direct observational data to 
guide us. Two basic approaches are available - full numerical simulations and semi-analytical/semi- 
numerical modeling. The simulations can handle the nonlinear dynamics, feedback effects and 
complex geometry, but are fairly computationally expensive. This drawback is being alleviated 
by innovative, efficient radiative transfer methods, some of which are able to efficiently utilise 
the latest Petascale computing facilities and specialist hardware (e.g. GPU accelerators). This 
now allows simulations which were impossible just a few years ago, with volumes and resolution 
matching both the intrinsic scales of the reionization process and the FOV and resolution expected 
for the SKA1 EoR experiment. Nonetheless, full simulations remain quite expensive and conse¬ 
quently a significant effort has been invested in the development of semi-numerical modelling and 
other faster modelling methods. These approaches are by nature more approximate, generally use 
simplified physics and are missing the non-linear effects. However, when they are carefully con¬ 
structed they still can retain many of the key features, while at the same time introducing the option 
of constructing large libraries of models to be used for interpreting the data, which is not feasible 
for full simulations. Much additional work is still required on these models, as well as on the nu¬ 
merical simulations, particularly in adding important physical processes which are still missing and 
in calibration and verification of semi-analytical models against numerical simulations. We have 
also outlined the requirements for SKA1-LOW-specific simulations and provided estimates on the 
resources which will be required for this effort. 
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